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Abstract 

We prove that the lubrication approximation is perturbed by a non-regular roughness of the 
boundary. We show how the flow may be accelerated using adequate rugosity profiles on the 
bottom. We explicit the possible effects of some abrupt changes in the profile. The limit system 
is mathematically justified through a variant of the notion of two-scale convergence. Finally, we 
present some numerical results, illustrating the limit system in the three-dimensional case. 



Introduction 

We study in this paper the effect of different very small domain irregularities on a thin film flow 
[ governed by the Stokes equations. There exist already some references on the subject. Roughly 
speaking, they may be ordered into three categories. The first one is devoted to the study of the 
roughness effects on the flow in a channel: the height of the channel, denoted h and depending on 
the horizontal component x, is fixed and the effect of small scales of the boundary is studied. In such 
■ studies, the height of the channel is written as (see [Hill]) 

I h{x) = /l2 (^x, . 

The second category is concerned with the specific study of the thin film assumption, that is when 
the height of channel is assumed to be small, of the form (see [9l [5l \T0\ lllj) 

i h{x) = ehi{x). 

The third category combines the mechanical point of view of the latter lubrication studies with the 
analysis of the roughness effects. Various limit models, in special regimes, are obtained depending 
on the ratio between the size of the rugosities and the mean height of the domain. In [3l HI HOl lllj - 
I the ratio is assumed to be of order one, namely 

I h{x) = eh2 (^x, — ^ 

and an asymptotic analysis is performed using an homogenization process. More recently, in [3], the 
authors study the case where the narrow gap is smaller than the roughness, namely 

h{x) = eh2 (^x, . 
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In all previous works, the ration between the two scales, modeled by the parameter a is such that 
a < 1. Moreover, the asymptotic model obtained in the case a < 1 is always the classical Reynolds 
equation. In [6J, the authors consider the particular case which does not enter in the previous 
framework: 

h(x) = eh2 {x, . 

They mathematically justisfy (through a variant of the notion of two-scale convergence) that an extra 
term modifies the standard Reynolds equation. 

There exists now a huge number of available data for rough surfaces, due to the increasing 
efficiency of the measurements (optic or laser technics, see e.g. [7]). On the one hand, the existence of 
multiple scales of rough oscillations is for instance emphasized in [13t [23] . But no analytical solution is 
available for complex geometries. These data are also often too cumbersome for numerical simulation. 
This motivates an approach using asymptotic analysis. On the other hand, the tribological literature 
does not provide any categoric answer for the choice of characteristic roughness parameters [151 121] . 
The aim of the present paper is thus to explore the effects due to abrupt changes in the rugosity profile. 
A direct application of this kind of study is the production of nano-scale electronic components which 
are nowadays formed by self assembly. In such a process, the solvent and the monomers are confined 
in a thin channel. The block copolymers create abrupt changes of the geometry (see for instance [12j). 

The case of a change of speed of order one could be easily deduced from our previous work |6J. We 
thus focus on more abrupt changes in the rough profile. In particular, we aim to compare the speed 
of the change, assumed of order e", a > 0, and those due to the characteristic speed of the rough 
oscillations, assumed of order e^, whereas e corresponds to the caracteristic height of the channel. 
We restrict ourself to the case where a < 2 to ensure that our limit model respects the structure of 
the classical Reynolds approximation [20]. 

We thus consider a model profile of the form 

h(x) =eh^ =e (/ii(2;) + e (^1 - (^)) (J)) ' < a < 2. 
Function -0 is introduced to mimic the change in the profile (see the Figure [T]). 




0(1) ©(e") 0(X) 

Figure 1: Domain Vl^ and the different scales 



The domain occupied by the fluid is then: 

= {(x, z) G E'^"^ xM, X £uj, < z < h{x)} 

where w is a domain of M'^^^. We choose d = 2 or d = 3 for the current applications. The change in 
the profile occurs in a subdomain Wj^ of cj. The domain uj is thus decomposed in u = cot U Wj- U uj^. 
We assume that the function ip belongs to L°°{uj) and is such that ip{x) = if x G and ^(x) = 1 
if X G Lv^. The positive function hi is the main order part of the roughness, while function /i2 
describes the oscillating part. For sake of simplicity, we assume that the fonction /il defined by 
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/iKx) = h2{x/£^) is an e^-periodic function. Actually, it would be sufficient for our needs to assume 
that /i2 is an admissible test function for the two-scale convergence (see Proposition 12. II below) . The 
first step of our method consists in introducing a new vertical variable. This change of variable 
motivates our last technical assumptions for the rough profile. We assume that ij) € C'^{W^~^), 
hi G H^{W^-^) and /i| G C^{T'^-^). 

Let us present our results. We prove that the change in the rough profile gives the same type of 
correction at the main order than the oscillating part /12, but only under an additional assumption 
linking the behavior of the roughness jumps with the one of the oscillating part of the profile. Without 
this additional assumption, the Reynolds approximation contains a non-explicit contribution (a kind 
of " strange term coming from nowhere" , see [8j ) . 

Theorem 0.1 Let uj he the horizontal projection of the domain of study. 

(i) Let us denote by ijp G L'^[u}]H^{T'^~^)) (respectively G L'^{uj) if a < 2) the two-scale limit 
of the jump approximation ^/>'^(x) = '4){x/e°') as e — )• and by Ub the velocity of the lower surface. 
The behavior of the thin flow is approximated by the following modified Reynolds equation on the 
pressure p of the fluid which only depend on the horizontal variable: 

div^ -^AV^p^ = div^ (/ii (SC/b - Q°)) on CO, 



,12 

where the functions A and B are defined on u by 



. 12 / ;v /2 /" N , \ 12 , ^ /2 Jn la ^ dtds 

N^^ Jo ) N^^ ' J^e^^^'/Us 

where N^{x) = j |Vx((l — 'i'^ {x , X))h2{X))'^ dX , j; G w. The "strange" function is non 
explicit. 

(a) Assuming moreover one of the following assumption, 

(HI) V{e'^xp^h2) strongly two-scale converges to V x{'^^h2), 

(H2) ip"^ strongly two-scale converges to tp^ , 
we recover a completely explicit Reynolds-type approximation: 



divJ^^I^V^pO) =dwJhiBUh 
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on UJ. 



Remark 0.1 R is crucial to clearly separate the oscillating profile /12 and the jumps characteristic 
function ip in such a study. Indeed, the limit behavior described in Theorem \0.1\ clearly depends on 
the link between the behavior of the roughness jumps with the one of the oscillating part of the profile 
(see Assumptions (HI) and (H2)). 

Remark 0.2 Let us describe a situation allowing the use of Assumptions (HI). Assume that ojj^ is 
a finite collection of disconnected subsets of u with measure of order e". As e 0, Uj, reduces to a 
finite discrete set of points {x^}i=i..Ar- Then 

N 
i=l 
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Since the Lebesgue measure does not see accidents, one easily checks with the definition of strong 
two-scale convergence in Proposition \2.1\ that 

N 
i=l 

This situation is illustrated through numerical exemples in Section^ 

Remark 0.3 We recall that the usual Reynolds approximation (corresponding to h2 = 0) reads 

div^(^^V^p°^ = diVa:{Y^b) on oj. 

The latter result is thus a perturbation of the classical Reynolds equation. The low order perturbations 
of the rough profile ( oscillations /i2 and abrupt change in the roughness ip ) both give a perturbation 
of the Couette component of the Reynolds model. This is a highly desirable effect for the applications 
of such a lubrication model. We explicitly recover a perturbation which was detected with formal 
computations by jlT^ , or more recently by 119]. 

Remark 0.4 The strong two-scale convergence mentioned in (ii) of Theorem \0.1\ could appear as a 
rather technical assumption. In fact, it means that the oscillation spectrum of the sequence belongs 
to the integer grid (see \22^ for a proof through Fourier analysis). The oscillation spectrum is then 
due only to the periodicity of the coefficients. 



The paper is organized as follows. In the first section, we give the model and the first estimates for 
the velocity and the pressure. In the second section we give the definition and the properties of a 
convenient two-scale convergence. The third section is devoted to the rigorous study of the upscaling 
process in view of proving Theorem 10.11 In Section [U we present some numerical results, illustrating 
the limit system in the three-dimensional lubrication context. 



1 Model and estimates 

1.1 Stokes model in lubricant context 

The vector field = {u'^,w'^) E x M (which describes the fluid velocity) and the pressure (given 
by the scalar function p^) satisfy the stationary Stokes equation^: 



-Au^-FV^/ = in^l^ (1) 

-A?i;" + a^/ = inf]", (2) 

diva;ii^ + d:,w^ = in 0% (3) 

= on dn^. (4) 



The viscosity of the fluid is set equal to 1 for sake of simplicity. 



*In all this document, the operators (like A) without index denote the operators with respect all the variables. To 
specify the operators with respect only one variable we use subscript notations, see for instance A^, or dz- 
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• In classical lubrication framework, the boundary conditions are the following: 

{{ub, 0) X ^ uj, z = 0, 

{euh{z),0) X G duj, 
(0,0) X euj, z = eh^, 

where Ub is a constant corresponding to the imposed velocity at the bottom of the mechanism, 
and Ub is a regular fonction (for instance an affine function) connecting Ub for z = and for 
z = . Notice that for such a problem (that is for an asymptotic study e — )• 0) the precise 
value for the velocity Ub is not primary. In fact, the physical quantity which persists to the 
limit is the total flux Ub{z) dz, see [3j. 

• Finally, the pressure is normalized by / = 0. 

We introduce a new vertical variable Z defined by Z = — to consider the same problem in the 

h 



fix domain 17 = w x (0, 1). Problem now read^ 



+ (t^ - ^)^^^^^ + 't^ • ^^^^^^^ - ^^'^^-^ - F^i-^ 



+VxP" - ^Zazp" = inJ7, (5) 
h 

+^dzp' = inn, (6) 

div^M^ - ^ • Zdzu^ + j^dzw^ = inQ, (7) 

= Ub on dn, (8) 

, A/i lV/i|2 ^. /V/i 

Observe that — —5— = div — — 

h h'^ \ h 

1.2 Estimates and dependance with respect to the small parameter e 

We begin by some estimates for the velocity. 

Lemma 1.1 There exists a constant C such that the velocity components satisfy the following uni- 
form estimates: 

l|V.C/"||(i2(o)). < J, (9) 

\\dzU\LHn-^y < C, (10) 

\\U'kLHn)r < C. (11) 

Proof: These estimates are directly derived from the original problem after an adequate lifting of 
the non-homogenenous boundary conditions. We write the standard energy estimate for Stokes type 
system and then use the change of variables to control the derivatives. We conclude with the Poincare 
inequality. □ 

We now give some uniform estimates for the pressure. 



^Rigorously, the function defined in the previous Stokes system ([T])-((l| is not egal to the function «^ used in this 
rescale domain. However, to not introduce numerous notations, we denote always the same. 
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Lemma 1.2 There exists a constant C such that the pressure satisfies 



I|V./||h-i(q) < (12) 

\\dzp'\\H-Hn) < f , (13) 
^^'lL^(n) < (14) 



Proof: Such estimates come from Equations ([S])-® estimated in H^^. This provides bounds on 
VxP'^ and dzp^ ■ Using Poincare-Wirtinger inequality (with the normalization hypothesis) we get the 
L^-estimates for p^. □ 

Notational convention: In what follows, for any function (f) € H^{ui x T'^"-'^), we denote by 
V0(x,x/e^) = Vcj)^ = Vx4>^ + e~^Vx</''^ its horizontal gradient. 



2 Definition and properties of a convenient two-scale convergence 

The proof of the homogenization process will be carried out by using a variant of the two-scale 
convergence introduced by G. Nguetseng in [18] and developed by G. Allaire in [2]. Let us give the 
basic definition and properties of this concept. 

Proposition 2.1 A sequence {v'^) of functions in L^(0) two-scale converges to a limit v^{x,Z,X) 
belonging to L^(f] x T'^^-'^), v^^v'^, if 

lim / v''{x,Z)^{x,Z,x/e^)dxdZ = [ [ v^{x, Z, X)'^{x, Z, X) dx dZ dX, 

for any test function ^{x,Z,X), X-periodic in the third variable, satisfying 

lim / \'^{x,Z,x/e^)\^ dxdZ = [ [ \^{x,Z,X)\^ dxdZ dX. 

Such a function ^ is called an admissible test function for the two-scale convergence. Note that Z is 
only a parameter for this definition. 

(i) From each bounded sequence {v^) in L^(ri) one can extract a subsequence which two-scale converges 

to some limit G L^(0 x T'^"^). The weak L"^ -limit of is v{x,Z) = fjd~i v^{x, Z, X)dX . 

(a) Let (u^) be a bounded sequence in L^(0, 1; H^{uj)) which converges weakly to v in L'^{0, 1; H^{uj)). 

Then v^-^v and there exists a function v'^ E L'^ {Q; [T'^^^)) such that, up to a subsequence, 

Vu^A Vt;(x) + \/xv^{x,y). 

(Hi) Let (v^) be a bounded sequence in L'^{Vl) such that {e^V ^v^) is bounded in {L'^{9)Y-^. Then, 
there exists a function € L"^ [0,; iT'^^^)) such that, up to a subsequence, 

v^^v^ and e^Vv^^Vxv^{x,Z,X). 

• Due to our assumptions for hi and /i2, functions hi{x), h2{x/e'^), dx^hi{x), 9a;. (e^/i2(x/e^)) = 
dx^h2{x / e^) can obviously be considered as admissible test functions for the two-scale conver- 
gence. 
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• We also note that the function V2 appearing in part (ii) of the previous proposition is the 
rigorous counterpart of the third term of the formal anisotropic expansion associated with the 
present setting: 

= Z, 4) + ev^{x, Z, 4) + e^^^(x, Z, 4) + • • • 

• An consequence of the above definition of the two-scale convergence is the following. 

Lemma 2.1 Let (u^) be a bounded sequence in L'^{Q,) such that (e^Wv^) is bounded in L?' , with rj < 2. 
Then the two-scale limit G L"^ [Q; [T'^^^)) of is such that 

Vxv° = 0. 

Proof: By (Hi) of Proposition 12. 1| e'^Vv^^S/ xv'^ . Since (e^'Vf^) is bounded in L^, it two-scale 
converges to some limit r]^ and (e^Vv^ = e'^~''^e'^Vv^) two-scale converges to 0. Thus V xv^ = 0. 

□ 

Note that the above lemma is a key result to treat the anisotropy of the rough profile. 

3 Rigorous derivation of the limit model 

3.1 Convergence results 

We infer from the estimates derived in Subsection II .21 the following result. 
Lemma 3.1 There exist limit functions 

such that 

e^p'^p^, u'^u°, Vxu° = 0, w'^w^ and Vxw^ = 0. 

Proof: In view of Lemma ll.ll we claim the existence of limits (up to a subsequence). Moreover, 
using Lemma |2.H we assert that S/xu^ = and Vxw^ = 0. □ 

Before passing to the limit in the equations, we state some auxiliary results. We begin by the pressure 
function. 

Lemma 3.2 The two-scale limit pressure is such that 

Vxp° = and dzp° = 0. 
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Proof: We multiply Equation ([5]) by an admissible test function in the form e^(j){x, Z,x/e'^) and we 

Ah iV/iP ,. fVh" 



integrate by parts. Unsing in particular — j-^ — = div ( -j— ) , we get 



Jn 

- [ e\V^hi+e~Hl-r)'^xhl-e^-''hlV^,r)-——-^———V^u'dz{Z^')dxdZ 
Jn hi + e{l - ip^Jh^^ 

- [ e\V^hi + e-\l-r)Vxhl-e^-''hlV^i^')- , ^ ^^,^ Zdzu'{V ^(f)' + \v x<^') dx dZ 
Jn hi + e{l - ip^jh^ 



+ / £ ieM.e\2 dzu' ■ dz[Z (p^) dx dZ 



/ 7T ^^TJJ2^zu' ■dz(l)'dxdZ - I e^p%£^d\w^(t)' + dbixct>')dxdZ 

+ I e\V^hi + e-\l-rWxhl-e^'^hlV^r)- , ^.^^ e^p'dz{Z(t>')dxdZ = Q. 

Jn hi + £[l - ip^jhl 

We bear in mind that the term is of lower order than if < a < 2 and of the same order if 
a = 2. Passing to the limit e — )• 0, we obtain 

lim / eVdivx^^dxdZ = / / p^dwxcj) dx dZ dX = Q. 
e^ojc JnJfi-^ 

Thus VxP^ = and the two-scale limit and the weak limit of (e'^p^) coincide. 
Now let G L'^{uj; H^{0, 1)). We write 

j^dz{£^p^)(l)dxdZ = x/fi ^0 

= - / e^p'dzcpdxdZ - [ p°dz(pdxdZ. 
Jn Jn 

It follows that dzp^ = 0. This ends the proof of the lemma. □ 

We now introduce auxiliary limit functions for sake of clearness in the computations. We define, 
7° G L2(J7;i7i(T'^-i)) and e° G L'^i^; H^{T'^-'^)) such that 

(1 _ ^e^usA ^2y^ ((i _ A y^^O, 

(1 _ V.^)VAeO, e2^,((^ _ ^.)2^.)A v^e°. 

Of course, if < a < 2, we have Vxt'^ = and VxS,^ = 0. We now study the two-scale limit of the 
vertical velocity component. 

Lemma 3.3 The vertical velocity component is such that 

Proof: We already mentioned that Vxw^ = (see lemma [3TT]) . We now prove that dzw° = 0. We 
multiply the divergence Equation d?]) by e(j){x, Z) where (/> S Integrating by parts, we obtain 

e 

t a 

n 



^ /ii + e(l - /ii + e(l - V^'j/il 

= 0. (15) 



hi +e(l - 
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Thanks to the relation 



H I r^TTT— 77urVx/i| • (1 - r)u'dz{Z<l>) = ! Vxh2i'^dx)dz{Z(P), 
£->o hi + e(l - Jq hi \Jjd-i ) 



which is justified by (1 — -0^)11^-^7'^ and to the relation 

.2-Q 



£^0Jn /ii + e(l - V^j/il 

^^ojf^ hi +e(l - '4)<')h\ 
+ l™n/ , ^ , .,^ /i|V^g^div.^^az(^<A) 



= / ^dz{Zcfy [ h2Vxl^dX, 
Jn "-1 Jfd-^ 

which is justified by e^VxH^ — ■iIj^)u^)-^ Vxt", we compute 

l™n(/" h J. ,,..e ^xh\-{\-^')u^dz{Z<fy- I ■ ^ " hl^x^l^' ■u'dz{Z<^) 

= j^^^dz{Z<p) I Vx{h2l'')dX 

= because of the periodicity of /i27*^- 



We thus infer from (|15|) that 

e^0jQhi+£{l-ljj^)hl J^hi 
It follows that dzw^ = 0. We conclude using the boundary conditions. □ 

Remark 3.1 We recall that the aim of this work is to get a perturbation of the Reynolds approxi- 
mation. And the previous result is characteristic of such a lubrication approximation: it states that 
we can neglect the vertical component of the velocity. 

3.2 Divergence equation 

Lemma 3.4 We have the relation 

/iidiv^n° - V^hi ■ Z dzu° + dziw^ +L) = 0, 

where vu^ = {1/d— 1) fjd-iWidX, with wi € L'^ {ft; {T'^^^)) is defined by the following two-scale 
convergence 

\ V{ew')^Vxwi, 
while dzL is defined by the following weak convergence 

dzL = hm dz ( Ze ^^ ~f''^'^ dwu') in V'{n) , 

e~^0 \ hi J 

L\z=o = L\z=i = 0. 
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Remark 3.2 The term is the rigorous counterpart of the first order term of a formal anisotropic 
expansion of w'^ : 

= w^{x,Z,x/e^) + ew^{x,Z,x/e^) + e'^w^{x,Z,x/e'^) H 

Proof: We multiply the divergence Equation ([7]) by a test function (f){x, Z) G We obtain: 

- [ u'-V,^+ [ -——^———V,hi-u'dziZ<P)dxdZ 
Jn J^hi + e{l-'il)^)hl 

+ / u ^ ,,,,e^^~'^'^ ^xhl-u'dz{ZcP)dxdZ 

- [ h ^ n le^ue ^Z^ y^r ■u'dz{Z(t))dxdZ - [ / -w'dz(l)dxdZ = 0. 

Jq hi + e(l - V^j/il e° 1 Jnhi + e{l - ^^)hl e 



It follows that 



Vx</'+ / ^Vxhi ■ ZdzvP(l)dxdZ 
Jn hi 



+ M I -^'~'^T'''}~'',ZT''^^ -^dzu^^d-dZ^ [ \ ' , ^^^.^dz^dxdZ 

= 0.(16) 

We note that l/(/ii +e(l — V'^)/!!) is a computable perturbation of l//ii for our convergences needs. 
Indeed, we have 



1 1 -e{l-ij^)k 



£ 



/ii +e(l - hi hi{hi + e{l - i;^)h^2)' 

On the one hand, using (1 — iIj^)'^u^^ and e'^Vxi^ — "0^)^"^^--^ VxC*^, we write 



/i2 Vxh2Zdz^°(t>dX 



lim 



+ ^™n/ /TTT^T^ — — --(i-^(x/e"))2v..(azn^)Z0 



hi 

f /,£2 2 

s-^oJq + e(l - V^j/il) 2 

= -^^ / ^xf'^^^zf^dX + J^^^ J hlVx{f)dz{Zcl)) dX 

= - J^^J^^ Vx/i2azC°^</>d^-^^y /i2 Vx{dzf)Z(l)dX 
= -J^^J Vx{hldze)ZcpdX 
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because of the periodicity of h2£,^. 

On the other hand, we define tu""^ by 



l™n/ h P^'e,he, ^"9z^= I (/ w'^'^- dX)dz<t> = - [ dz ( [ w'^^^dx)^, (17) 



which means that (1 - ^l}^)h\w^ lh\ w'^'^^'dX weakly in L^{n). Then Relation reads: 



/ M°-V^0+ / -^V:,./ii • Z9zuV + [ dz( [ w''''''dx)(P = -A- B (18) 



where 



We write 



f 1 I 

B = lim / -—w^dz4'- 

e^o Jq e hi 



A = lim [ ^V,(e(l - i;')hl) ■ {Zdzu')(t> 
hm(/ £i^^div..^9z(Z0) - / <'-f^^h .hi.u^-dziZ^) 

+ /^^^f^n^-V.(a.(Z0))) 

= lin, / £(l^mdiv.^^az(z<A) 

Jq hi 

^. /■e(l-V')^^. ^ /• e(l - V')/i2^. 

^^^Jn hi ^^oJq hi 

= - lim / dz (Z^ii^^div.n^) + lim / ^i^^div...^^. (19) 
^^oja ^ hi ) e^CijQ hi 



The first term in the right hand side will give the conservative contribution dzL to the final divergence 
equation. Let us compute the last term. We multiply the divergence Equation d?]) by the test function 
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e{l-ij'')hl(j){x,Z). We get 
lim / , ^, w Oz0 

= lim f ed\w^u'{l - r)hl(i) + lim / ^iLzif)/|^v,/ii • u'dz{Zct^) 

e^Q e~^0 hi + e(l - 

+ lim f ^i_^Vx/i| • (1 - i^'fu'dziZcj)) 

- lim / ^^^^V^r ■ (1 - r)u'dz{Zc^) 
= lim [ ediv,n^(l-V")/i|0+ [ ^ [ V x{hl) ■ ^ dziZcj)) dX 

r 2 

-^™n/ , ^,(1 -{l-r?diy.u^dz{Zct^) 
'^^oJq hi +e{l -tp^)h^2 2 

= \\r^J^ediY,u'{l-il^')hl4> + j^^J Vx{hl)-edz{Z<P)dX + j^^J hi V dz{Z<t>) dX 

= lim / ediv,n^(l-V")/i|0+i [ ^ [ V xihl^) dz{Z(p) dX 
e-^o 2 J^hi J 

= lim / e div^n^ (1 — ip^) ho (j) 

because of the periodicity of h2£,^. 

Bearing in mind (|17p . we infer from the latter relation that: 



e^oJ^hi{hi+e{l-i^^)h'2) ^^oj^hi ^ ^ 

dz( [ ■u^""'^dx)</. (20) 
\Jjd-i J 



because of the definition of w""^^ . 

f 1 1 

Let us now compute B = lim / ——w'^dz4'- We write 



e^oj^ £ hi 

1 1 



-^w^dz4> 
n^hi 



f 1 1 

Jn ^ hi 

^ l^div{I{x/e^))iew'^^dz^) 
^ I{x/e^) ■ V{ew'^dz(t>) 



d-l Jn hi 
where I is the 1-periodic function such that I{x) = x if x € T"'^^. Then 

B = --\[([ X-Vxwidx)^dz4 
d-l Jn^JT'i'i ^ hi 



^ ^ widx)-^dz(t> (21) 



d- Jn ^7x^-1 ^ hi 

f — dz(t) - - f ^^^^ 
Jn hi ^ Jn hi 
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where = fjd-i wi dX, the function wi € L'^{Q; H^{T'^ ^)) being defined by 

Using (|18|) -(|21|). we obtain the announced result. □ 

The conservative form of the Hmit divergence equation obtained in the previous lemma is sufficient 
to get an explicit Reynolds approximation of the flow, even if function L is not explicitly computed 
(see [6]). By the way, we end this section by computing L under an additional assumption. 

Lemma 3.5 Assuming moreover one of the following assumption, 
(HI) V('(/'^/i2) strongly two-scale converges to Vxifp^h2), 
(H2) ip^ strongly two-scale converges to ip^, 
(H3) a<2, 

then L = and the limit divergence equation reads 

/iidiv^.n° - V^/ii • Zdzu^ + dzw^ = 0. 
Proof: From the assumptions (HI) or (H2), it follows that 

h2{Xh\x,X) = /i2(X)(l - ip\x,X))uo{x,Z), 

h2{X)e{x,X) = h2{X){l - V'°(x,X))'no(x,Z). 

The two scale limit of Z e(l - ^/;^) /i| divn^ is hiL = Z (1 - /i2 div^u^ where € L2(J1; iji(T'^-i)) 
is the anisotropic two-scale limit defined by 

Let us compute divxu^- To this aim, we multiply the divergence Equation ([7]) by e(/){x, Z,x/e'^) 
where (p is an admissible test function for the two scale convergence. Integrating by parts, we obtain 

/ edivu'cP'- [ [{l-i^')Vxhl-e^-'^hlV^il^')-Zdzu \ f 
Jn Jn /ii + e(l - 

/ii + e(l - 7n /ii + e(l - 

We recall that edzu^^ (a strong two-scale convergence due to the fact that is bounded in 

{L^{Sl)f) and w^^Q (see Lemma l3.3p . We also have 

hm / ((1 - i,')Vxh\ - e2-"/j|V,V') • Zdzu \ ^ ,f ^ = / / -3- divx(/i2^5z7°) </> d^, 

because 

((1 - V^)Vx/i| - e^~"/i|V^VO^i^ = e^div^.((l - V'^)/i|n^) -e^{\- V'^)/i|div^'u^ 

^ divx (/i27°) ■ 

We thus conclude that 



divxni = -^divx (/i2 • Zdzl""). (22) 
hi 
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With the additional assumption (HI) or (H2), we infer from (1221) and from V^jti" = that 

divxu^ = ^divx(/i2(l - V'°)^9zn0) = -ivx(^2(l - /)) • Zdzu'^. 
hi hi 

We conclude by computing 

/ LdX = ^dz{Z^dzuo)( [ /i2(l-V'°)Vx(/i2(l-V'°))^i^) =0. 

Jfd-i h^ ^Jja-i J 

Assume now (H3) . Following the lines of the proof of Lemma 13.31 one easily check that the 
assumption a < 2 leads to 

(1 - V^)w;^^0. 

Now, using the test function e(\ — V'^)(/)(x, Z, x/e^) instead of x/e^) in the latter derivation 

of divxti^, we state that 

2 1 „ .o^ 



e(l - V^)div'u=^ 7-divx(/i2 • Zdz^ 
hi 



Integrating over T , we recover once again L = 0. □ 
3.3 Momentum equation 

Lemma 3.6 Letq^ be the "strange" function defined by the following weak convergence in L'^{uj-, H^"^ (0,1)): 

|V(e2(l - l^Zdzin' - n') + Zdz[zdz{\e'p'))v{e\l - V'^)/i||2) - 

The limit of the first momentum equation is 

N^Zdzu"^ - + ^?V,pO + g° = 0, 

where 

i(x)= / \Vx{{l-i^\x,X))h2iX))\^dX, xeuj. 



Remark 3.3 Assuming moreover (HI) or (H2), that is V(^/'^/i2) strongly two-scale converges to 
Vx{tp^h2) or ip^ strongly two-scale converges to , and using dzp^ = 0, one checks that q^ = 0. 

Proof: Let us pass to the limit e ^ in the momentum equation. We multiply Equation ([5]) by 
e^(/)(x, Z) where (p € Hq{uj; Hq{0, 1)) and we obtain: 

Jn hi + e{l - il)^)hl 

Jn (/ii + e(l - 



/q (/ii + e(l - 



"-1 + e(l - ^ )^2 



14 



Passing to the limit e — ?> 0, we get 



lim 

e^O In V hi 



e-i-0 



hi 



0. 



(23) 



We express the first term of the left hand side of (f23]l using once again ([5]), but for the test function 



,,ii-r)hi 



e dz{Z4>). For this choice of test function, we obtain: 

hi 



-3 



= / , ^ n ,e,ue {'^-hi+e-\il-r)'^xhl-e^-'-h',V,i;'))-Zdzu'v( ^^ P^h z{Zcl>) 
Jnhi + e{l-ip^)h^2 ^ hi 

+ / , , f ,,..e iy-hi + e-\l-rWxhl-e^-''hlV^r)-eV,u'dz(z ^^~P^'^ dz{Z<t^) 
Jnhi+e{l -'ip^)h^2 ^ hi 



ik.^e(l-Ah^^ --' eOzu^ . 8z[Z^^^MZ,) 



Passing to the hmit e — > in the latter relation, we get: 



lim / V.n--£^V( ^^ P^' dz{Zct>)) 
e^Oj^ \ hi I 

= lim I ^^{(\-^-)^^h\-e^-'^h\V^r)-Zdzu'vC'^^^^dz{Z^)\ 
^^^jQ.hi^ \ hi J 

+ lim / e3pMivf^i^^^9z(Z0)) 

e^^J^ \ hi ) 

- hm / ^eV((l - ^')^xh\ - e2-"/i|V.V0 • dz[z^^—^^dz[Z4>)\. (24) 

e-^Oj^hi ' \ hi J 
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Inserting the previous computation in (j23|l . we obtain: 



■lim / -^((1 - V")Vx/^| -g^-^/^IVV'") • Zdzu'V(e^ ^^ f ^''h z{Z(t>)) 
"-1 ^ hi 

-lim [ e'p'div( ^^~P'^h z{Z<t>)) 



hi 

+ lim /" e3/div(^i^^^^5z(Z</.)) = 0. 

We thus note that the combined use of ([23]) and (pl|) . that is the use of the momentum equation ([5]) 
for two different ad hoc test functions, is sufficient to pass to the limit without studying the limit 
behavior of the non-bounded sequence (ep^). This is a subsequent improvement of [6j. At this step, 
we have proven that: 

-lim [ ^\{l-^^)Vxh%-e^-''h%Vi)^?Zdzu^dz{Zc^) 

+ lim [ J-£V((1 - i;^)Vxh^ - e^-^h%vr) ~ P^h z{ZdziZ^)) 

+ hm / ^ ^ ' \^ 'i—L^dzu^-dz{Z^4)) 

+ / T29zu^ ■Oz'p- [ P°dw^<P = 0, 

that is, combining the first and third lines of the latter relation, 

lim [ ^e^p-{{l-^')Vxh%-e''-"hlvr)^^^—^^dz{Zdz{Z(t^)) 

e^Oj^hi hi 

Jn hf 

+ / 13^^"° ■ ^^'^ - / P°div,0 = 0. (25) 

It remains to exhibit the partial differential equation corresponding to the latter weak formulation. 
Lemma 13.61 is proven. □ 

Remark 3.4 Assuming a < 2, the definition of the "strange" function becomes: 

I V(e2(l - pzaz(n^ - u°) q° weakly in L'^{uj; H-'^{0, 1)). 

Indeed, if a < 2, lim(eVV(e2|(l-V'=)/i||2)) = lim(eV(l-V'^)^V(e2|/i||2)) = /^,_, T]°Vx{\h2\'^)dX 
where rp is the two-scale limit of the sequence (e^p^(l — ip'^)'^). This limit is computable without 
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assuming the strong L'^ -convergence of (e^p^) as in A first proof is in 116^ - Here we simply use 
the two-scale convergence. With the slight modification of the proof of Lemma \3.Sl we obtain 

lim ( I £2/(1 _ divx0^ dxdZ- I e^p^e^'^'V^' (p" dx dz) = 0, 

for any admissible two-scale test function cp^ . Since a < 2, it follows that V xff = 0. Thus, by the 
periodicity of /i2, 

I rfVx{\h2?)dX = rP [ V x{\h2?)dX = {). 
Jfd-i Jja-i 

3.4 A modified Reynolds approximation: proof of Theorem 10.11 

Let us check that the system formed by the equations presented in Lemmas 13.41 and 13.61 corresponds 
to a "modified Reynolds" system. The hmit momentum equation is (see lemma [32]) 

N^Zdzu^ - dlzu^ + hlV^p^ + g° = 0. 
We integrate it in the variable Z (for each fixed x) and find, with the boundary conditions: 

u\x,Z) = ( [' /%^'^(-)(^^-*^)/2(V./(x) + P^)dtds 
\Jo Jo hi{xy 



Ub. (26) 



Next, integrating with respect to Z the divergence equation (see lemma [3 

div^.(/iin°) + dz{wi +L- ZV^hi ■ m°) = 0, 

and taking into account that the velocity li^i + L — ZV^/ii • cancels for Z = and for Z = 1 we 
obtain 



div^ hiu^dZ^ = 0. 



With the previous expression (|26|) for the velocity the following pressure equation is obtained 

div, (/i?AV,p°) = div, {hi{BUi, - Q°)) 
where A and B are two functions defined in w by 
12 



A{x) 



Ln^{x)/2 / -N^{x)ty2 ,U_A_ (Mx)/2 _ . Jo Jo ^ ^ ^ 

N^{x)\! io J ^ J^' e^A^)^y^ ds ' 

B(x) = - (e^^^^'^/'^ - 1) - X G u 

assuming that function y i— > (e^/^ — l)/y is continuously extended to and Q*^ is defined by 

QO(a;)= /Y r f e''^^-)^''-'')l\\x,t)dtds 
Jo \Jo Jo 



1 

JO 
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Part (i) oi Theorem 10. II is proven. Part (ii) is straightforward since the hmit momentum equation is 
now given by: 

N^Zdzu'^ - dlzu"" + /i?V,p° = 0. 
We thus have g*^ = and then = 0. 

4 Some numerical illustrations 

In this section, we present some numerical results. The goal is to quantify the effect of roughness on 
a portion of a surface. So we compare the solutions of the classical Reynolds equation (corresponding 
to the smooth domain) and those obtained from the equation with partial roughness, see Theorem lO.il 
on page [3] and Remark 10.21 after. 

Physically, we are in the framework of lubrication. Consequently, we are interestd in the pressure 
forces in a flow between two surfaces whose relative velocity is imposed. As usual, all the geometry 
of the gap between the two surfaces is taking into account by the upper surface (it corresponds to 
the height h\) and the relative velocity is taking account by the lower surface (its imposed velocity is 
denoted by C/fe, while the upper surface is at rest). Moreover, we know that the problem is well posed 
if and only if we give boundary conditions: either by imposing Dirichlet condition on the pressure, 
either by imposing Neumann type flow. A physically interesting situation is the following: we give a 
Neumann condition at the "entrance" of the channel (that is a flux denoted by Qe), and we impose 
a pressure at the other boundaries. 

More precisly, for the next simulations we use the following data (see also Fig. [2]) : 

a; = [0,1] X [0,1], /ii(a;,y) = (2a; - 1)2 + 0.5, f7fe = (l,0) and Qe = 0.5 



4.1 Without rugosities 



To situate the next results, the first figure that we present (see Fig. [21 on right) corresponds to 
the case of flow in a smooth domain. In that case, we solve the classical Reynolds equation, see the 
Remark 10. 3 1 We use the Freefem++ language (see http:/ /www.freefem.org/ff++/ ) and consequently 
a variational form of the Reynolds equation: for all test function q vanishing on the boudaries {x = 1}, 
{y = 0} and {y = 1}, we get 



/j3 



Qeq = 0. 



{x=0} 



Finally, we use the Scilab software (http://www.scilab.org/ ) to view results. 



4.2 Some rugosities 

From the results of Theorem lO.il it seems simple to simulate the rugosities effets on the upper surface. 
Morally, it is sufficient to specify first: on which area of the surface roughness ago, and secondly: 
quantify the roughness using their form. These data allow to evaluate the function N^, and therefore 
the coefficients A and B involved in the statement of the Theorem 10.11 For standard rugosities, 
for instance described by a function of the form acos{f3X), it is reasonable to take the following 
examplarity values 

N^{x,y) = 0, A{x,y) = l and i? = 0.5 on smooth part, 
N^{x,y) = 2, A{x,y) = 1.08696 and B = 0.58739 on rough part. 
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Bottom velocity 



Figure 2: On the left figure, we draw the three-dimensional domain and indicate the main data. On 
the right figure, we plot the pressure profile without rugosity, that is the solution to the classical 
Reynolds equation. 

We then present three cases where the domain a; is partially rough (see Fig. [3l H] and [5]) . We can 
note two remarkable points. First, the effect induced by the roughness is not a local effect. Thus, 
although roughness occupy only a portion of the domain, then the pressure is disrupted throughout 
the domain, including upstream roughness. Secondly, the effect of roughness is quite significant even 
if these wrinkles are present only a small part of the domain (see Fig. [5]). 




Figure 3: On the left figure, we indicate the part of the top boundary which is rough. On the right 
figure, we compare the pressure profile without rugosity (grid) and the pressure profile obtain with 
the corresponding rugosities (plain). 

Acknowledgments - The second author of this work has been partially supported by the ANR 
project n° ANR-08-JCJC-0104-01 : RUGO (Analyse et calcul des effets de rugosites sur les ecoulements) . 
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Smooth part Rugosity part 

Figure 4: On the left figure, we indicate the part of the top boundary which is rough. On the right 
figure, we compare the pressure profile without rugosity (grid) and the pressure profile obtain with 
the corresponding rugosities (plain). 
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^ Smooth part Smooth part 

Figure 5: On the left figure, we indicate the part of the top boundary which is rough. On the right 
figure, we compare the pressure profile without rugosity (grid) and the pressure profile obtain with 
the corresponding rugosities (plain). 
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